Universal pulse sequence to minimize spin dephasing in the central spin decoherence 

problem 
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We present a remarkable finding that a recently discovered [G. S. Uhrig, Phys. Rev. Lett. 98, 
100504 (2007)] series of pulse sequences, designed to optimally restore coherence to a qubit in the 
spin-boson model of decoherence, is in fact completely model-independent and generically valid for 
arbitrary dephasing Hamiltonians given sufficiently short delay times between pulses. The series 
maximizes qubit fidelity versus number of applied pulses for sufficiently short delay times because 
the series, with each additional pulse, cancels successive orders of a time expansion for the fidelity 
decay. The "magical" universality of this property, which was not appreciated earlier, requires that 
a linearly growing set of "unknowns" (the delay times) must simultaneously satisfy an exponentially 
growing set of nonlinear equations that involve arbitrary dephasing Hamiltonian operators. 

PACS numbers: 03.67.Pp; 03.65.Yz; 76.60.Lz; 03.67.Lx 
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Introduction. The spin-boson (SB) model, where the 
environment is modelled by a bosonic bath of simple har- 
monic oscillator quanta, is a famous and extensively used 
general technique for studying the quantum decoherence 
of a system coupled to an environment [l[ . The model is 
often used in the context of quantum dissipation and de- 
coherence analyses of wide classes of couplings between a 
system and its environment. Recently, Uhrig 0] discov- 
ered a series of 7r-pulse sequences that optimally, with 
respect to the number of pulses, decouple a qubit from a 
bosonic bath in the SB model, thus protecting the qubit 
from decoherence. Because errors are inherently intro- 
duced by each pulse, the promise of optimal dynamical 
decoupling (DD) is of great significance to the intense 
recent activity in quantum information processing where 
minimizing qubit decoherence and successfully carrying 
out quantum error correction protocols are crucial. 

A particularly important issue in view of the remark- 
able power of the Uhrig DD (UDD) in fighting qubit de- 
coherence is its applicability in realistic situations beyond 
the idealized SB model context of its discovery ^] . To un- 
derstand the general applicability of these sequences, we 
set out to apply it to a drastically different model of de- 
coherence, the central spin decoherence problem; specif- 
ically, we consider the spectral diffusion (SD) process, 
quantum dephasing of an electron spin qubit coupled 
to a slowly fluctuating nuclear spin environment (a spin 
bath). Our investigations led us to discover that UDD 
transcends all models and is an optimal decoupling se- 
quence for any dephasing Hamiltonian when delay times 
are sufficiently short because, with each additional pulse, 
it kills successive orders of a time or Magnus [1, [j| ex- 
pansion 0. The apphcability of the UDD sequence to 
the SD model, with its extreme contrast to SB, already 
indicates its powerful generality. What we find is much 



more general - we show that the Uhrig sequence is model- 
independent. 

How to preserve the state of a qubit in a bath is an 
important theoretical and practical consideration in the 
field of quantum information. A large energy splitting 
between the qubit's two logical states (e.g., through the 
application of a magnetic field for a spin qubit) com- 
pared with the temperature of the bath may result in 
long relaxation (Ti) times; however, the relative phase 
of a superposition state may not be preserved by this 
strategy so that dephasing {T2) decoherence ensues. As 
a strategy to combat dephasing and an example of DD, a 
sequence of 7r-pulses may be applied in order to rapidly, 
on the time scale of the system dynamics, flip the qubit 
in between time intervals of free evolution (this may be 
generalized for any quantum system using inverting Q 
pulses as a generalization of 7r-pulses). In the simplest 
case, the Hahn spin echo occurs after applying a single 
TT pulse midway through the system's evolution. Con- 
catenated fvl DD (CDD) sequences can successively im- 
prove coherence times, but at the considerable overhead 
expense of exponentially increasing the number of ap- 
plied pulses. UDD performs its magic, not only in the 
SB model but whenever delay times are sufficiently short, 
with a mere linear scaling in the number of applied pulses. 

Pulse sequence echoes. Considering only dephasing de- 
coherence, the effective Hamiltonian (for any model) may 
be written in the form TL = ^^|±)H±(±| with |+) and 
— ) as the two qubit ket states. For a given pulse se- 
quence with intervals Ti between successive pulses, the 
evolution operator is then U = ^_|_|±)C/±(±| (or U — 

'^±\-f)U±{±\ if there are an odd number of pulses) with 



exp I ~iHipT2 ) exp I —iH±Ti 



(1) 
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In order to characterize the coherence decay, we con- 
sider the transverse component of the qubit's expecta- 
tion value; normalized to a maximum of one, the pulse 
sequence echo, ve, is defined in this way such that 
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{W) where W = ULU+, the (...) 

denotes an appropriately weighted average over the bath 
states (we use equal weights justified for temperatures 
large compared to nuclear Zeeman energies), and ||...|| is 
the magnitude of the resulting complex number. 

Uhrig series. The UDD sequence with n pulses may 
be defined by [2| 



Aj - 1) 



t, (2) 



ioi 1 < j < n + 1 [corresponding to tj in Eq. ([T|)] where 
t is the total sequence time. This series was shown 0] 
to optimally, with respect to number of pulses, suppress 
decoherence in the SB model. We find, remarkably, that 
for any form of dephasing Hamiltonian, with no assump- 
tions about 7i±, this sequence yields ve = l — O (t^"+^) . 
Equivalently stated n pulses in the UDD series re- 
moves the first n orders of the Magnus expansion of u^;! 
In comparison with CDD, the number of required pulses 
scales exponentially with respect to the orders of the 
Magnus expansion that are cancelled (It is impor- 
tant to note, however, that CDD can compete with UDD 
when the Magnus expansion does not converge well) . We 
will return to a discussion of the universality (Hamilto- 
nian/model independence) of the UDD series at the end 
of this Letter after we demonstrate its consequence for a 
model that is drastically different from SB. 

Spin boson versus spin bath. In the SB model, the spin 
qubit interacts with a bath of non-interacting bosons: 
H = J^i'^M^i + + ^i)' where b represents 

boson operators and Sz is the z spin operator for the 
central spin. In contrast, the spin bath model treats in- 
teractions of a central spin, such as a localized electron 
in a solid, with interacting nuclear spins in a solid-state 
lattice. Exemplifying the spin-bath model, SD is a de- 
phasing of the central spin as a result of fiuctuations of 
the bath-induced effective magnetic field caused by intra- 
bath interactions. In the limit of a large applied field, the 
polarization of the central electron spin and the nuclear 
spins must be individually preserved (the electron having 
a gyromagnetic ratio that is typically 2000 times larger 
than those of the nuclei), so that the Hamiltonian is 
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(3) 

Spectral diffusion: cluster expansion. Despite the 
mesoscopic size of the solid-state baths that typically 
contribute to SD, often involving many millions of nu- 



clear spins, it is feasible to compute ve — 



{W) 



with 



problem into manageable sub-problems that each involve 
only a few nuclei. (This expansion was successfully ap- 
plied [§| to the problem of Si:P donor electron SD yielding 
remarkable agreement [9| with experiments 0.) Con- 
sider expanding W such that W = X]n=o where 
iyt"l contains contributions to W that involve n separate 
clusters of "operatively" interacting nuclei. To be spe- 
cific, the set of nuclei involved in a term of W^^^ must all 
be connected together via factors of bilinear interaction 
operators to form a single connected cluster. Clusters 
have spatial proximity when interactions are local. If it is 
possible to approximate (M^^^l) by only including clusters 
up to some small size that is much less than the number 
of nuclei in the bath, N, and if the initial bath state 
is effectively uncorrelated (e.g., a random bath), then 
(iy["l) « {W^^^y^/n\. In this "cluster approximation,' 



VE 



(W) 



exp ^Re |(iyt^l)|y This approximation 

is extremely useful because it is possible to treat (lyl^l) 
perturbatively in cases where the perturbation would fail 
for (W) directly due to the vast numbers of multiple clus- 
ters involved in {W) (by definition W^^^ involves only sin- 
gle clusters). In our calculations, we approximate (lyl^l) 
using a perturbation in cluster size (we include contribu- 
tions from successively increasing clusters until conver- 
gence is achieved). The justification for this perturbation 
in cluster size is that each additional nucleus in the clus- 
ter requires an additional bilinear interaction factor. An 
expansion in cluster size rather than a direct expansion 
in orders of intra-bath interaction factors (e.g., diagra- 
matically [ll| ) is simply more convenient. 

Spectral diffusion: time perturbation. In addition to 
the intra-bath perturbation, it is also possible for (ly'^l) 
to converge in a time expansion. The time perturba- 
tion applies when the tj are small compared to all in- 
teractions time scales of the system. The time pertur- 
bation is only really relevant if it is applicable on the 
time scale of the decay (e.g., T2 observed in spin echo). 
Different clusters operate, in the sense of contributing to 
SD decay, on very different time scales depending largely 
upon differences in the HF interactions (A„) among clus- 
ter nuclei, and these interactions are inhomogeneous over 
the bath. When there are enough clusters with enough 
influence operating on the shortest time scales so that 
these clusters dominate the decay, the decay time may 
be small compared to all interaction time scales such 
that the time perturbation is relevant. This, in turn, 
depends on the distribution of the An determined by the 
shape of the electron wavefunction. With these consider- 
ations and assuming that intra-bath interactions are lo- 
cal (e.g., dipolar interactions), electrons in quantum dots 
with gaussian-shaped wavefunctions will tend to, in gen- 
eral, exhibit short-time behavior SD on the decay time 
scale while donor-bound electrons, with radially expo- 
nential wave functions, will not l^]. As an example of 



W = ULU+ using a cluster expansion that breaks up the a situation in which the time perturbation is not appro- 
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FIG. 1: (color online). Various UDD (left) and CDD (right) 
pulse sequences in side-by-side comparison where red lines 
mark the time of each applied n pulse. The sequences are 
numbered by their order which corresponds, respectively, to 
the number of pulses (UDD) or the level of concatenation 
(CDD). Orders of UDD and CDD with the same number of 
pulses are lined up horizontally. 



priate, the Hahn echo decay of donor-bound electrons in 
Si:P has the form exp(— r^ "^) Q which cannot be ex- 
plained by any Magnus or time expansion. The Hahn 
echo decay of GaAs quantum dots, on the other hand, 
exhibit exp (— r**) behavior Q explained by the lowest 
order of a time perturbation expansion. 

Spectral diffusion: pulse sequences. Using the cluster 
expansion technique, we are able to test the effective- 
ness of DD strategies that use 7r-pulse sequences in the 
real-world SD problem. We have previously theoretically 
verified [l3| the powerful effect of SD suppression when 
applying concatenations of the Hahn echo sequence de- 
fined recursively by p; :— Pi_i — > tt — > pj_]^ — > tt with 
Po T. These CDD sequences cancel successive per- 
turbative orders of both the intra-bath [l^l and Mag- 
nus/time 0, 0, expansions with each concatenation. 
The main advantage of this sequence compared with 
UDD is that it operates on the intra-bath perturbation 
which can be applicable on a much longer time scale, for 
the time between pulses, than a time expansion (intra- 
bath coupling is on the order of ms while HF interac- 
tions limits the time expansion on the order of /is) . Each 
concatenation, however, essentially doubles the number 
of applied pulses, leading to exponential overhead. The 
main advantange of the UDD series [Eq. ^] is that it 
yields successive time expansion cancellations with each 
added pulse for the SB model, a linear overhead! Fig- 
ure [1] shows a side- by-side comparison of UDD and CDD 
sequences. 

Spectral diffusion: results. We show cluster expan- 
sion results (to the lowest non-trivial order in the intra- 
bath perturbation), with coupling constants of Eq. ^ 
obtained using models described in Ref. Q . We compare 
the effects of CDD and UDD sequences on the coherence 
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FIG. 2: (color online). Numerical results showing the nat- 
ural logarithm of echoes from CDD and UDD sequences for 
Si:P donor electrons in natural Si. Red (green) denotes UDD 
(CDD) sequence results; the Hahn echo decay is shown as 
a solid black curve for reference. Sequences with the same 
number of pulses have the same line pattern. UDD has no 
particular advantage (in the high fidelity regime) to CDD in 
Si:P, but UDD with an even number of pulses performs better 
than those with an odd number. 



of both an electron bound to a P donor in natural Si 
[Fig. [2] and a quantum dot electron in GaAs [Fig. [3], 
plotted as a function of the total time t of one iteration 
of the pulse sequence. In Si:P, where only the intra-bath 
perturbation (and not the time perturbation) is applica- 
ble, CDD maintains high fidelity (e.g., 10"'' decay) for a 
longer time than UDD. The CDD sequences cancel out 
perturbative orders of the intra-bath coupling with each 
concatenation fl3j; however, UDD sequences cancel only 
the first and, if there are an even number of pulses, sec- 
ond orders. In the GaAs quantum dot, in which time per- 
turbation is valid, CDD still does better (maintains high 
fidelity longer) than UDD for equal cancellations of order 
[Fig. [3] (a)]. However, for equal numbers of pulses, UDD 
preserves coherence far better than CDD [Fig. [3] (b)]. 
In this case, UDD cancels out perturbative orders of a 
time expansion (exhibited by the successively increasing 
slopes of the curves) with each extra pulse, exactly what 
it was shown to do in the SB model. The only signifi- 
cant difference between Si:P and the GaAs quantum dot 
is the applicability of time perturbation expansion due 
to their respective electron wavefunction shapes. Despite 
the stark, qualitative difference between SB and SD mod- 
els, UDD proves to be optimal, in the short time limit 
(where accessible), for both. 

Universality. What is most striking, however, is that 
UDD is generically optimal for cancelling orders of a time 
perturbation for any dephasing decoherence! With full 
generality, let Ti.± = Xq ± Xi without making any as- 
sumption about the commutation properties of Xq and 
Xi. This Hamiltonian, along with the pulse time inter- 
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FIG. 3; (color online). Numerical results showing the natural 
logarithm of echoes from CDD and UDD sequences for quan- 
tum dots in GaAs with a 25 nm Fock-Darwin radius, 8.5 nm 
quantum well thickness. Red (green) denotes UDD (CDD) 
sequence results; the Hahn echo decay is shown as a solid 
black curve for reference, (a) Sequences of the same order, 
having the same initial slope, use the same line pattern, (b) 
Sequences with the same number of pulses have the same line 
pattern. Note that the Hahn (black) and CDD (green) curves 
in (a) and (b) are the same, shown for comparison with UDD 
curves. 



generically true for all integer n, and numerical calcu- 
lations (up to n = 14) are consistent with this assertion. 
This finding is remarkable! The UDD series satisfy an ex- 
ponentially growing set of C... = non-linear equations 
with degrees of freedom (jj) that merely grow linearly 
(one for each pulse). 

Conclusion. We conclude by emphasizing our key find- 
ing that the UDD [§| sequences restore coherence op- 
timally and generically in a model independent man- 
ner through successive cancellations in orders of a 
time/ Magnus Q perturbation expansion with a low over- 
head of a single pulse for each order of cancellation. 
While we assume ideal, instantaneous pulses in our anal- 
ysis, careful pulse shaping can theoretically[l5[ mitigi- 
tate effects of the finite width in experimentally applied 
pulses. The universal Hamiltonian-independent applica- 
bility of this series of pulse sequences, originally proposed 
within the narrow constraint of a spin boson model, is 
simply miraculous. 

We thank Professor Gotz Uhrig for bringing his op- 
timal sequences with linear scaling (in the number of 
pulses) to our attention, and Lukasz Cywihski for ex- 
tremely helpful advice in the presentation of this paper. 
We thank Professor Deepak Dhar for bringing Ref. [u\ to 
our attention. This work is supported by LPS-NSA and 
ARO-DTO. 



vals, Tj, for a given pulse sequence will determine the U± 

operators that go into w^; {W) with = TJ~ U^. 

Taking the real part of W yields a more convenient ex- 
pression, Re |(W^)| = 1 - (AtA)/2 with A = U+ -U-. 

Since > Re it only serves to make a more pes- 

simistic estimate to ignore the imaginary part of (W); we 
may therefore restrict our attention to the real part in our 
perturbative analysis. 

Expanding the exponentials of U± [Eq. ([T])] and collect- 
ing terms with common sequences of Xq^i operators, we 
define Ci^^i^^,,, as the coefficient of a term in A = U^ — IJ^ 
for a corresponding sequence of operators such that 



A 



OC 1 

E E 

m— 21,22,...,*^.— 



^Xi^Xi^...Xi, 



(4) 



These C... coefficients may be written as a sum of prod- 
ucts of the Tj pulse delay times, defined by Eq. ^ for a 
sequence of n pulses, by expanding Eq. ([T|). It so happens 
that for any m < n, all Cij^ij^...^^^ ~ O (t™) coefficients 
are identically zero! We have explicitly proven this for 
(up to) n < 9 using computer integer arithmetic by repre- 
senting each C... as a polynomial of a = exp {m /(n + 1)), 



exploiting the fact that a 



n+l 



-1, and not- 



ing that two or more points placed symmetrically (i.e., 
equally spaced) around a circle in the complex plane sum 
identically to zero. We conjecture that this result is 
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